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Abstract. In this paper we present an algorithm, called conauto-2.0, that can 
efficiently compute a set of generators of the automorphism group of a graph, 
and test whether two graphs are isomorphic, finding an isomorphism if they are. 
This algorithm uses the basic individualization/refinement technique, and is an 
improved version of the algorithm conaiito, which has been shown to be very 
fast for random graphs and several families of hard graphs [8]. In this paper, it is 
proved that, under some circumstances, it is not only possible to prune the search 
space (using already found generators of the automorphism group), but also to in- 
fer new generators without the need of explicitly finding an automorphism of the 
graph. This result is especially suited for graphs with regularly connected com- 
ponents, and can be applied in any isomorphism testing and canonical labeling 
algorithm (that use the individualization/refinement technique) to significantly 
improve its performance. Additionally, a dynamic target cell selection function is 
used to adapt to different graphs. The resulting algorithm preserves all the nice 
features of conauto, but reduces the time for testing graphs with regularly con- 
nected components and other hard graph families. We run extensive experiments, 
which show that the most popular algorithms (namely, nauty [9,10], bliss [6,5], 
Traces [12], and saucy [2,7]) are slower than conauto-2.0, among others, for the 
graph families based on components. 



1 Introduction 

The Graph Isomorphism problem (GI) has been puzzling computer scientists for several 
decades, because it could not be shown to have polynomial complexity, but it could not 
be proven to be NP-complete either. However, testing if two graphs are isomorphic is 
useful in a number of contexts, ranging from chemistry [3,16] to computer vision [1]. 
Hence, the design of algorithms and the construction of tools that are able to solve the 
problem efficiently for a large collection of problem instances has significant interest. 

Related Work For the last three decades, nauty [9,10] has been the most widely used 
tool for graph isomorphism testing and canonical labeling. However, Miyazaki proved 
[11] that nauty required exponential time for a family of colored graphs. McKay noted 
also that nauty would also require exponential time for unions of strongly regular 
graphs. All this has encouraged researchers to develop new tools, to try to overcome 
this drawback. Some of them are based, like nauty, on canonical labeling. Examples 
are bliss [6,5], Traces [12], and nishe [15,14]. Another tool, named saucy [2,7] solves a 



related problem: computing the automorphism group of a graph. It is specially designed 
to efficiently process big sparse graphs. 

A different way to tackle the GI problem was suggested in [8]. The tool developed 
in that work, called conauto, does not generate a canonical labeling of the graphs being 
tested, but instead looks for similar sequences of vertex partitions. To do so, it uses a 
limited search for automorphisms on the graphs. This algorithm has good performance 
in practice for different families of graphs, like for example the graphs of Miyazaki. 
However, its inability of computing the whole automorphism group restricts the benefit 
obtained from known automorphisms. That reduced its performance with some families 
of graphs, like Latin square graphs or the point-line graphs of Desarguesian projective 
planes. One major contribution of the foUowup work [13] was a way to avoid back- 
tracking when processing union graphs, which could be applied in a more general field. 
Some ways to improve algorithm conauto were also suggested as open problems; in 
particular, computing the full automorphism group of the graphs, and recording and 
using information on non-isomorphisms to prune the search space^. 

Contributions In this paper we describe a new algorithm, called conauto-2.0 that, using 
the basic approaches of [8], extends the functionality of the algorithm, and significantly 
improves its performance. The main contribution is a new theorem that is applicable 
under very simple conditions, and prunes the search for automorphisms. This has a 
large impact in many classes of graphs, like, for instance, graphs built from connected 
components. This result is directly applicable to any other algorithm that uses the indi- 
vidualization/refinement approach to compute the automorphism group or the canonical 
labeling of a graph. 

One key issue of every automorphism group computation or canonical labeling al- 
gorithm is cell selection. In many cases, the target cell chosen for individualization 
conditions the ability of the algorithm to effectively find the generators of the automor- 
phism group. This is particularly critical in the case of Desarguesian projective planes. 
In fact, knowing how an algorithm selects the target cell enables the construction of 
graphs that make an algorithm exponential in time. That was the approach of Miyazaki 
against nauty. While nauty and bliss use a static cell selector. Traces uses a dynamic cell 
selector that uses multi-refinements. However, although this approach shows especially 
powerful for non-Desarguesian projective planes, it is too expensive for simpler graphs. 
In conauto-2.0, a dynamic cell selector is used but, for efficiency, it is not isomorphism 
invariant and, hence, cannot be used for canonical labeUng. However, it is able to choose 
the best target cell in the case of many different graph families (for example the Desar- 
guesian projective planes, or Miyazaki' s graphs). Finally, two additional features of 
conauto-2.0 are that, in addition to test for isomorphism, it computes tlie complete 
automorphism group of a graph (set of generators, size of the automorphism group, 
and orbits), and it uses detected non-automorphisms in a similar way to the failures 
management of bliss [5]. 

^ The latest versions of bliss (versions 0.65 and above) [5] apply similar ideas: component recur- 
sion, and use of failures to prune the search space. However, bliss-0.72 stops with an internal 
error when processing some graphs (based on unions of connected components). The authors 
of bliss have been notified of the problem. 



We have carried out extensive experiments to compare the practical performance 
of conauto-2.0 with nauty-2.4, bliss-0.35, and bliss-0.72'* to test for isomorphism with 
several graph famiUes. (Only a small subset of results has been included here.) Our 
experiments show that conauto-2.0 is the fastest for most graph families, and when it 
is not, its performance is similar to that of the others. We have also evaluated the per- 
formance of several programs to compute the automorphism group of several selected 
graphs. The results show that the best performance is presented by blissO.72, Traces, 
and conauto-2.0, without a clear winner among them. 

Structure Next section defines the basic notion of sequence of partitions on which 
conauto-2.0 is based. Then the key features of the algorithm are described. In Section 4, 
the main theorem is presented. Finally, a practical performance evaluation is done. 



2 Definitions and Notation 



A directed graph G is the pair {V, R) where F is a finite non-empty set of vertices 
and i? is a binary relation. The elements of R are called arcs. An arc {u, v) G R is 
oriented from utov. We use the term graph to refer to a directed graph^. For the sake 
of clarity we may not be fully formal in the following definitions. More details and 
formal definitions can be found in [8,13]. 

An isomorphism of graphs G = {Vg, Rg) and H = {Vh, Rh) is a bijection be- 
tween the vertex sets of then, / : Vg — > Vh, such that (v, u) G Rg iff {.f{v), ,f{u)) € 
Rh- Graphs G and H are isomorphic, written G ~ H, if there is some isomorphism of 
them. An automorphism of G is an isomorphism of G and itself. 

Given a graph G = {V, R), we assume R to be given as an adjacency matrix 
Adj{G) = M with size |F| x |y| as follows. 



if {u, v) ^ RA {v, u)(^R 1 if (m, v) (jt R A {v,u) € R 
2 if {u, v) € RA {v, u)iR 3 if {u, v) G R A {v, u) e R. 



Given a vertex v £ V and a subset V C V, the available degree of v with V is the 
3-tuple (£>3, D2,Di) counting the number of vertices in V connected to v with the 3 
types of adjacencies 3, 2, and 1, respectively. If all the elements of the tuple are zero 
then, V is disconnected from V . The available degree is used in our algorithms to order 
the vertices of G. 



Sequences of partitions The conauto [8,13] algorithms all use the same basic approach 
to test for isomorphism between graphs G and H: they build a sequence of partitions 
for one of the graphs and try to find another sequence of partitions for the other graph 
with the same underlying structure (compatible). A sequence of partitions defines an 
ordering of the vertices of the graph. Then, if compatible sequences of partitions are 
found, the corresponding orderings yield the isomorphism between G and H. 

* Two versions of bliss have been considered because bliss-0.72, although faster than bUss-0.35 

in some cases, sometimes crashes, as was mentioned. 
' An undirected graph is a graph whose arc set R is symmetrical, i.e. (u, w) £ J? iff (w, m) £ -R. 



Let us consider a graph G = {V,R). A. partition of a set 5 C 1/ is a sequence 
<S = {S\, Sr) of disjoint nonempty subsets of S such that S = lJi=i ^i- The sets Si 
are called the cells of <S. A sequence of partitions of G is a list of partitions <S°, <S*. 
The sequence starts with the trivial partition = {V), and partition 5*^+^ is obtained 
from partition using some refinement. Each refinement can be a set refinement or a 
vertex refinement (or vertex individuaUzation). In both cases a pivot cell is chosen. In a 
set refinement, each cell of the initial partition 5*^ is divided into smaller cells in S'^^^ 
as a function of the available degree of each of its vertices with respect to the pivot 
cell. In the vertex refinement, a vertex p of the pivot cell is chosen, and each cell of the 
initial partition is divided into smaller cells as a function of the adjacency of each of its 
vertices with respect to p. Each refinement can discard some vertices from the partition, 
either because they have been used as pivot vertex in a vertex refinement, or because 
they have no more Unks with the rest of vertices in the partition. The sequence ends 
when all remaining cells are singleton. In a sequence of partitions, level k refers to «S*^ 
and its associated parameters (e.g., type of refinement, and pivot cell used). 

Consider two graphs G and H, and sequences of partitions >S° , . . . , <S* and T" , . . . , T* 
of them. We say that the sequences are compatible if both have the same number of lev- 
els t, at each level the corresponding partitions are compatible as well (i.e., have the 
same number of cells, pairwise with the same size and the same available degree be- 
tween them), the refinements applied at each level k is the same, and the pivot cell is in 
the same position, and the final partitions <S* and T* have the same adjacencies between 
cells. The following theorem shows that having compatible sequences of partitions is 
equivalent to being isomorphic. 

Theorem 1 ([8]). Two graphs G and H are isomorphic iff there are two compatible 
sequences of partitions <S°, <S* and T°, for graphs G and H respectively. 

The challenge to derive fast isomorphism testing algorithms based on sequence of 
partitions is to deal with backtracking in the search for a compatible sequence of par- 
titions. Backtracking occurs when vertex refinement is used in the original sequence, 
and the pivot cell has more than one vertex. Our algorithms always guarantee that this 
happens only when the partition has no singleton cells and it is not possible to refine 
the partition by means of a set refinement. When this happens at level k, we say that we 
have a backtracking point at that level, or that fc is a backtracking level. 

In order to deal with backtracking points, the conauto algorithms use vertex equiv- 
alence and automorphism detection. In a graph G = {V, R), two vertices u,v G V are 
equivalent at the level I of a sequence of partitions if there is another sequence of parti- 
tions, compatible with the first one, in which the first I partitions are the same. Observe 
that two vertices that are equivalent at level / are equivalent at all levels k,k < I. Two 
vertices equivalent at some level belong to the same orbit. 

3 Description of the Algorithm conauto-2.0 

Algorithm conauto-2.0 compares two graphs G and H for isomorphism. It works in 
the following way. First it computes a sequence of partitions for graph G. Then, unlike 
previous versions of conauto which only perform a partial (polynomial time) search 



for automorphisms, conauto-2.0 performs a full search for automorphisms. This yields 
a complete set of generators for the automorphism group of the graph. The automor- 
phisms found allow removing backtracking points from the sequence of partitions (a 
backtracking point at a level I is removed, when all the vertices of the pivot cell in 
level I are equivalent at that level). If, after removing backtracking points, the result- 
ing sequence of partitions of G still has backtracking points, a sequence of partitions 
for graph H is generated, and the search for automorphisms and backtracking point re- 
moval in that graph is performed. Note that, if the sequence of partitions of G does not 
have backtracking points, then it is not necessary to generate the sequence of partitions 
forF. 

Then, a backtracking process is used to find a match between the vertices of the two 
graphs. In this process, the sequence of partitions of one of the graphs is used as target, 
and an attempt is made to find a compatible sequence of partitions for the other graph. 
This uses the knowledge of the automorphism groups of the graphs, what drastically 
prunes the search space. The sequence of partitions with less backtracking points is 
chosen as the target for the match. Observe that if the target sequence of partitions has 
no backtracking points, the matching process has no backtracking. 

Pivot cell selection. When the sequence of partitions of a graph is computed, at each 
backtracking point, a pivot cell must be chosen for vertex individualization. In conauto- 
2.0, at a level k, for each possible combination of cell size and available degree, the first 
cell in the partition, with that size and available degree, is considered. The first vertex 
in each of these cells is individualized and the obtained partition is subsequently refined 
until an equitable partition is reached*". For each vertex x, assume the equitable partition 
is reached at level l^. If the partition at that level is a sub-partition of the partition at 
level k, then the cell of x is chosen without testing any other. Otherwise, we sum the 
number of discarded vertices and the number of cells of this resulting partition. Then, 
the cell of the vertex x which yields the biggest such sum is chosen as the pivot cell. 
Note that, since only one vertex in each cell is considered, the choice is not isomorphism 
invariant and, hence, cannot be used for canonical labeling. However, it works very well 
in practice for computing the automorphism group of a graph. 

Recorded failures. In [13] one of the authors of this work suggested that recording 
failures during the search for automorphisms could help pruning the search space. Re- 
cently, the same idea has been proposed by Junttila and Kaski in [5]. The use of failures 
in conauto (version 2.0) and bliss (version 0.65 and above) are similar, but in our case 
they are used both in the search for automorphisms and during the matching process. 

Search for automorphisms. The search for automorphisms is performed iteratively 
traversing all the backtracking points of the sequence of partitions from the last to the 
first. At each backtracking point, every vertex of the pivot cell (except the pivot vertex 
p used in the original sequence of partitions), which has not been found yet to be equiv- 
alent to p, is chosen for vertex individualization. It is used to generate an alternative 

* Partition S = {Si,...,Sr) is equitable if, for all i,j £ {!,..., r} and all m, v £ Si, the 
available degree of u and v with Sj is the same. 



sequence of partitions using the same pivot cells and refinement procedures used in the 
original sequence of partitions. If a compatible sequence of partitions is found, the cor- 
responding automorphism is recorded as a generator, and the vertex equivalences are 
updated accordingly. If all the vertices of the pivot cell at some level are found to be 
equivalent, then the backtracking point is removed. 

The search for a compatible sequence of partitions that induces an automorphism 
of the graph is performed by a backtracking process that explores every feasible path in 
the search tree (all feasible pivot vertices at each backtracking point traversed). Known 
automorphisms and sub-partitions (as described below) are used to prune this search. 

Use of sub-partitions to prune the search for automorphisms. Consider the following 
definition. 

Definition 1. Let k and I, k < I, be two backtracking levels. We say that is a sub- 
partition ofS'^ if h such that i ^ j, Sj C Sf^, and 5j C S^. (I.e., each cell ofS^ 
is included in a different cell ofS'^.) 

As it will be proved in next section, when the search for automorphisms is being per- 
formed at level k, it is enough to consider partitions up to level I, where I is the first 
level reached such that 5' is a sub-partition of 5*^. If a compatible alternative sequence 
of partitions can be generated up to level I, then it is possible to infer an automor- 
phism without the need to generate any further partition. Observe that the property that 
partition 5' is a sub-partition of partition S'^ can be evaluated before the search for 
automorphisms starts. 

4 Sub-partitions Theorem and Correctness 

In this section we present the main result that has been used to improve the performance 
of conauto-2.0. It is expressed as Theorem 2 below. 

Let S = (5°, ...,S*) be a sequence of partitions for graph G = {V,R), where 
<S' = {S\,...,Sl.) and denotes the the set of vertices in partition S^, for all i e 
{0, t}. We consider two backtracking levels k and I so that >S' is a sub-partition of 
.S*^ (see Definition 1). Let p be the pivot vertex used for the vertex refinement at level 
k. Consider a vertex p{of the same cell) that, if used instead of p as pivot vertex at 
level k, generates an alternative sequence of partitions T'^+^j .... T' that is compatible 
with S'''^^ , . . . , »S' . We use to denote the set of vertices in partition . 

We define now the following sets. Let E = V'' \V' md E' = V'' \ be the 
vertices discarded in the original and alternative sequences of partitions from level k to 
level I, respectively. Then, let us define A = EdE' bc\hc common vertices discarded, 
B = E\Athe vertices discarded only in the original (sub)sequence, C = E' \Athe 
vertices discarded only in the alternative (sub)sequence, and D = V^O the common 
vertices not discarded. Clearly, E = Au B, and E' = Au C. Let us also, for each 
X e {E, E', A, B, C, D}, define the subset Xi = X n S^. Observe that \Ei\ = \E'i\, 
and hence \Bi\ = \Ci\. 

As it was observed in [13], the graph induced by E is isomorphic to the graph in- 
duced by E', and there is an isomorphism of them that matches the vertices in Ei to 



Algorithm 1 Map the vertices of V with vertices of W 

ComputeSubPartition{G, S, T, V\ W^) : void 

1 for each Vi G F ' do Wi Vi end for 

2 for all f4 G C do 

3 J ^ I 

4 while ej G j4 do j <— A; : Cj = e'^. end while 

5 Wi ej where I : wi = el 

6 end for 



those in i?-, for all i. Furthermore, let ei, 62, e\E\ and e'^, ejg| be the vertices 
in E and respectively, in the order in which they have been removed from their 
respective sequences of partitions. Then, mapping Ci to gives such isomorphism. The 
(sub)sequence of partitions S\...,S* also imposes an order among the vertices of VK 
Let ui, W2, be such order. Algorithm 1 defines an order wi, W2, .... w^y^ of the 
vertices of so that the mapping of Vi to Wi, for each Vi G V\ extends the isomor- 
phism induced by the sequences of partitions to an automorphism of graph G. What 
Algorithm 1 does to construct the order wi, W2, W|yi| is to start as vi, V2, 
and replace each vertex in this sequence that belongs to C by a different vertex from B. 

Theorem 2. Let T' and be two subpartitions ofS'^, compatible between them, and 
k < I, the identity mapping for vertices in V \ V'^, combined with the mapping between 
E and E' induced by the sequence of partitions, and the mapping between and 
obtained with Algorithm 1, define an automorphism of graph G. 

5 Performance Evaluation 

In this section we compare the performance of conauto-2.0 against other algorithms 
for graph isomorphism testing and automorphism group computation. The experiments 
have been carried out in an Intel i7 Q 720 @ 1 .6GHz with 8GiB of RAM under Ubuntu 
10.04. All programs have been compiled with gcc 4.4.3 with their respective default 
configuration, but modified to perform isomorphism testings or automorphism group 
computation, depending on the experiment. The first experiment considers some graph 
families, and is intended to see how the size of the graphs affect the running time of 
isomorphism testing programs. Each point shown in the plots corresponds to the average 
running time of 100 executions with different instances of the corresponding graph. 
In the second experiment, some singular cases have been considered, to compare the 
running times for automorphism group computation. The times shown are for a single 
instance of each graph. The CPU time Umit of each execution has been set to 10,000 
seconds for both experiments. 

Many graph families are used by the different graph isomorphism and canonical 
labeling algorithms developers. Unfortunately, there is not enough room here for all 
of them, and not all of them are significant. In Figure 1 we show the results we have 
obtained with a few selected families. They are of different types, dense and sparse, 
they are all hard for most algorithms. The families considered are the following. 
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Fig. 1. Performance evaluation for testing isomorpliism. 



- Kronecker Eye Flip graphs (KEF). They are part of the benchmark of bliss [4]. 

- Cubic Hypo-Hamiltonian graphs clique-connected {CCH_cc). This family is pro- 
posed by the authors. It is built using as basic block two non-isomorphic Hypo- 
Hamiltonian graphs with 22 vertices. Both graphs have four orbits of sizes: one, 
three, six, and twelve. A graph CHH_cc{m, n) has n complex commponents built 



from m basic components. The components of a complex component are connected 
through a complete m-partite graph using the graphs that belong to the orbits of size 
three of each basic component. The n complex components are intercormected with 
a complete n-partite graph using the vertices of each complex component that be- 
long to the orbits of size one in the basic components. With this family, bliss-0.72 
returns an internal error for many benchmark graphs of 88 or more vertices. 

- Latin square graphs. These are ^3(71) latin squares of order n, for non-prime n. 
They were hard for the previous versions of conauto, but not for conauto-2.0. 

- Miyazaki Augmented2 graphs . These are taken from the benchmark of bliss [4] . For 
this family of graphs, bliss-0.72 gives a significant improvement over nauty-2.4 and 
bliss-0.35, but its performace is far from that of conauto-2.0 for large instances. 

- Point-line graphs of Desarguesian projective planes. This family was extremely 
hard for previous versions of conauto. However, thanks to the adaptive selection of 
the pivot cells, it has now a performance comparable with both versions of bUss. 

- Tripartite graphs, partially connected. These graphs are proposed by the authors, 
and are built from two 13 vertices directed tripartite graphs. Four vertices of each 
component have an arc to three vertices of each other component. Thus, each com- 
ponent is connected to all the other components of the graph. For graphs of 78 
vertices or more, bhss-0.72 generates internal errors in many cases. There are not 
results for nauty because for some instances on 26 vertices, nauty-2.4 is not able to 
finish within the time limit. 

- Unions of strongly regular graphs. The components are strongly regular graphs of 
29 vertices each. Each vertex in one component is connected to every vertex of the 
other components. Thus, this family is extremely dense. Observe that only conauto- 
2.0 is able to finish within the time limit for graphs of more than 203 vertices. 

- Paley toumaments. For this family of graphs all algorithms exhibit nice time com- 
plexity. However, conauto-2.0 is clearly the fastest of the four programs tested. 

To conclude, we can say that conauto-2.0 outperforms all the other algorithms or has a 

similar behavior for the graph families considered. 

Figure 2 shows the results of the second experiment. As mentioned, in this ex- 
periment we test the performance of computing the automorphism group of different 
graphs. All the graphs used have been taken from the benchmark of bliss [4]. As can be 
seen, the best performance is presented by blissO.72, Traces, and conauto-2.0, without 
a clear winner among them. 
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execution has exceeded the time hmit. 
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A Proof of Theorem 2 



The proof of following lemmas can be found in [13, chapter 7, p. 70-75]. 
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Fig. 3. Partition of into subsets Ai, Bi, Ci, and A. for all i. 



Lemma 1. Let M = Adj{G). For each u S E,foralli S {1, ...,ri},for all v,w S 5*-, 
Muy = Muw and M^„ = M^u- Similarly, For each u' G E', for alii G {1, r;}, for 
all v', w' e TI Mu'v' = Mu'w' and My'^' = Aft«'u'- 

Lemma 2. For each i,j G {1, .... r;}, there is some a € {0, 3} such that for all 
u G Bi, V £ Ci, w e Di, u' e Bj, v' e Cj, and w' G Dj, M^v' = M^w' = Myu' = 

Myyjl = M^u' = Myjy, = aUCl Mu'v = My^lyj = My^y = M y> y, = Myjly = 

My,,y = a-^ (We define 0"^ = 0, 1"^ = 2, 2"! = 1, and 3-^ = 3.) 
Let us define two families of partitions of Ai for i,j G {1, r}: 

Af = {xGAi:yuGBi,v'G Cj, M^y, = Myy,} 

= {X G Ai -.MU G Bi, V' G Cj, M^y, ^ Myy,} 

Note that, since the vertices of A^ are unable to distinguish among the vertices of Cj, 
then, if M^v' 7^ Myy' for some u G Bi or some v' G Cj, then M^v' 7^ Myy' for 
all u G Bi and all v' G Cj. Hence, each pair of sets A'^-' and A"-* defines a partition 
of Ai. Note also that, since each vertex in Ai has the same type of adjacency with all 

the vertices in U U Di (from Lemma 1), then for all x G A'^-' , u G Bi, v G Ci, 
w G Di, u' G Bj, v' G Cj, and w' G Dj, M^y, = M^y, = M^y,, = Myy' = M„^/ = 
Myyi = Myyji = Myjyi = Myjyi (f rom Lcnmia 2) . 
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Fig. 4. Partition of Ai into subsets A^, and Af. 



Lemma 3. For all i € {1, r}, let = fj^^i A':\ and let Af = IJ^^i A"\ Then, 
any isomorphism o/Ge andGs' that maps Ge^ to Ge'., maps the vertices in Af among 
themselves. 

Lemma 4. Gb is isomorphic to Gc, and there is an isomorphism of them that matches 

the vertices in Bi to those in Ci, for all i G {1, ...,r}. 

Lemma 5. Gyi and Gyyi are isomorphic, and there is an isomorphism of them that 
maps the vertices in Sj to the vertices ofT^forall i G {1, r;}. 

Let I be a backtracking level and let 5' = {S{, 5^) be the partition at that level, 
then observe that, for all i G {1, r}, Ggi is regular. 

Once we reach level I, generating an alternative sequence of partitions, which is 
subpartition of a level k < Z, we apply Algorithm 1 due to obtain an automorphism of 
G. 

Observation 1 A vertex discarded up to level I has the same adjacency with all the 
vertices of each cell of level I. Otherwise it would have split such cell. 

Observation 2 The vertices that belong to the same cell at level k and have not been 
discarded up to level I, must belong to the same cell at level I, since the partition of 
level I is a subpartition of the partition of level k. 

Observation 3 For all i G {1, |iJ|}, Cj and belong to the same cell in level k. 
This follows from the compatibility of the sequences of partitions. 

Observation 4 Algorithm 1 substitutes a vertex that belongs to a cell in level k with a 
vertex that belonged to the same cell at that level k. This follows from Observation 3. 

Proof of Theorem 2: 

Let xi, x\y\yk\ be the order induced by the sequences of partitions on the 

vertices discarded up to level k. Then x\,..., x\y\vk\,e\, ...,e\E\,vx, ...,v\v'^\ is the or- 
der induced by the sequence of partitions S. Let us call this order base. Analogously, 
a;i, .... .X|y\yi |, e'^, .... e'^^^.wi, .... w^y'i is the new generator computed. We now prove 
that they define an automorphism of the graphs, i.e.: 

1. Foreachi,jG{l,...,|y\l/'=|},M,.,^. = M^,^, . (Trivial). 

2. For each i,j G {1, Me.e, = M^,^y 

3. For each i,j e {1,...,|F'|},M„.^^ =M^,^.. 

4. For each i e {1, \V \ V'']},,] G {1, \E\}, M,,e, = M,^,,. 

5. For each i G {1, \V \ V'lj^j G {1, \V'\}, M,,,^ = M,,^.. 

6. Foreach^G{l,...,|S|},.?G{l,...,|^/'|}.A^,.., = M,>^^. 

Case 2 follows directly from the compatibility of the sequences of partitions. 

Case 4 follows directly from the compatibility of the sequences of partitions. 

Case 5 follows from the fact that, if no discarded vertex of a sequence of partitions 
has spUt a cell, that is because that vertex had the same adjacency with all the vertices 
of this cell. Since the partitions at level / are subpartitions of the partition at level k 



and Observation 4, all the vertices discarded before level k have the same adjacencies 
with all the vertices of each cell of level I. Otherwise, those cells would have been split 
before level k. 

From Observation 4, the vertices of a cell of level I that belonged to come cell of 
level k are substituted by Algorithm 1 by vertices that also belonged to that cell of 
level k. Hence, in level I they belonged to corresponding cells. Hence, from this fact. 
Lemma 1 and the compatibility of the sequences of partitions, Case 6 follows. 

From Lemma 5, there is an isomorphism that maps the vertices in corresponding 
cells of level Z. It is easy to see that mapping the vertices of D to themselves, we obtain 
one automorphism of Gd (which maps the vertices of Di among themselves for all i). 
We only need to prove now that having substituted the vertices of Ci by those of Bi 
preserves the adjacencies among them, and from them to the vertices of D. 

From Lemma 2, the vertices of Bi have with the vertices of Dj the same adjacen- 
cies that have the vertices of Cj with those of Dj, for alH, j. Hence, the substitution 
preserves these adjacencies. 

Finally we need to prove that Algorithm 1 preserves the adjacencies among the 
vertices of B and the vertices of C. 

From Lemma 3, we know that the vertices of B and C are not mapped to any vertex 
of A"'. Consider first the vertices of B that are mapped by the isomorphism induced by 
the sequences of partitions to vertices of C. Each of these vertices of B have the same 
type of adjacency with the other vertices of B than the corresponding vertex of C with 
the other vertices of C. 

Let us consider now the vertices of B and C that are mapped to vertices of A (in 
fact, from Lemma 3, we know that they are mapped to vertices of A^). In this case. 
Algorithm 1 maps each vertex bi G Bi to a vertex a e Ci such that the isomorphism 
induced by the sequences of partitions maps bi to some vertex Ui in A'^, and this vertex 
Gi is mapped to Cj (there could be a chain of vertices in ^4? instead of only one, but that 
makes no difference in the argument). Since the vertices of A? have the same type of 
adjacency with all the vertices of Bj, Cj, for all j (from the definition of A^). Hence, if 
vertex bi was mapped to a^, then it had the same type of adjacency with all the vertices 
of Bj, Since was mapped to c,, then Cj also had the same type of adjacency with 
all the vertices of Cj, (the same that Ui had). Hence, vertex Ci had the same type of 
adjacency with all the vertices of Cj, than vertex bi had with all the vertices of Bj, for 
all j. Thus we prove that the mapping generated by Algorithm 1 yields an isomorphism 
of Gb and Gc, completing the proof. 



